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Abstract 

The task of calculating operator dimensions in the planar limit of M = 4 super Yang- Mills 
theory can be vastly simplified by mapping the dilatation generator to the Hamiltonian of an 
integrable spin chain. The Bethe ansatz has been used in this context to compute the spectra 
of spin chains associated with various sectors of the theory which are known to decouple in 
the planar (large- A^ c ) limit. These techniques are powerful at leading order in perturbation 
theory but become increasingly complicated beyond one loop in the 't Hooft parameter A = 
gy M N c , where spin chains typically acquire long-range (non-nearest-neighbor) interactions. 
In certain sectors of the theory, moreover, higher- loop Bethe ansatze do not even exist. We 
develop a virial expansion of the spin chain Hamiltonian as an alternative to the Bethe 
ansatz methodology, a method which simplifies the computation of dimensions of multi- 
impurity operators at higher loops in A. We use these methods to extract previously reported 
numerical gauge theory predictions near the BMN limit for comparison with corresponding 
results on the string theory side of the AdS/CFT correspondence. For completeness, we 
compare our virial results with predictions that can be derived from current Bethe ansatz 
technology. 
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1 Introduction 



A two-dimensional parameter space emerges around the BMN limit of = 4 SU(N C ) super 
Yang-Mills (SYM) theory and the dual pp-wave limit of IIB superstring theory on AdS§ x S 5 . 
This space is parameterized on one axis by the perturbative gauge theory expansion in 
powers of the 't Hooft coupling A = gy M N c , and on the other by a string theory curvature 
expansion away from the Penrose limit in inverse powers of the string angular momentum 
J. The original comparison made by Berenstein, Maldacena and Nastase [Tj, which has 
sparked a number of direct tests of the AdS/CFT correspondence in recent years, lies at 
the intersection of the one- loop order (0(A)) gauge theory prediction and the zeroth-order 
(0(J°)) limit in the string curvature expansion. To explore a larger region of this duality 
landscape, a number of studies have pushed the string theory calculation to higher orders in 
the 1/ J curvature expansion j2j EJ E] , while operator dimensions in the gauge theory have 
been computed to higher orders in A (see, eg., EH 13 El E] ) • 

The first 1 / J curvature correction to the fully quantized string theory near the pp-wave 
limit of AdS 5 x S 5 was studied for two-impurity string states in El H|- In this setting, 
the extended superconformal multiplet of the theory is a 256-dimensional multiplet built 
on a scalar primary (or "highest- weight" ) state. This analysis was extended to the three- 
impurity, 4,096-dimensional supermultiplet in ^H]- To test predictions of the AdS/CFT 
correspondence, it was necessary in the course of the latter study to obtain higher-impurity, 
higher- loop operator dimensions from the gauge theory near the BMN limit. This is hard 
to do by standard diagrammatic methods, but the problem has been drastically simplified 
by the recent discovery that, in certain sectors of the gauge theory, the dilatation operator 
can be mapped to the Hamiltonian of an integrable spin-chain system. Calculating operator 
dimensions is therefore equivalent to finding the eigenvalue spectrum of these spin chains 
and certain established techniques associated with integrable systems, most notably the 
Bethe ansatz, have proved useful in this context (for a general review of the Bethe ansatz 
method, see ^1]). The utility of this approach in the setting of M = 4 super Yang-Mills 
theory was first demonstrated by Minahan and Zarembo ^2]- For operators with two R- 
charge impurities, the spin chain spectra can be computed exactly via the Bethe ansatz. For 
three or higher- impurity operators, however, the Bethe equations have (to the best of our 
knowledge) only been solved perturbatively near the limit of infinite chain length ^21 E31 El- 
Furthermore, at higher- loop order in A, the spin chain Hamiltonians typically acquire long- 
range or non-nearest-neighbor interactions for which a Bethe ansatz may not be available. 
For example, while the action of the spin chain Hamiltonian in the "closed su(2|3)" sector 
is known to three-loop order [7j, the corresponding long-range Bethe ansatz is not known 
(though it may well exist). A long-range Bethe ansatz does exist for the particularly simple 
"closed su(2)" sector of the theory [HIE], and our methods will provide a useful cross-check 
on these approaches to higher-order gauge theory anomalous dimensions. 

To improve on the current limitations of Bethe ansatz techniques, we have developed a 
virial approach to the spin chain systems of Af = 4 super Yang-Mills theory. The generic 
spin-chain Hamiltonian acts on single-impurity pseudoparticles as a lattice Laplacian and 
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higher iV-body interactions among pseudoparticles are suppressed relative to the one-body 
pseudoparticle energy by inverse powers of the lattice length L. Surprisingly, this expansion 
of the spin chain Hamiltonian is truncated at 0(L~ 3 ) in certain subsectors of the theory, 
allowing straightforward eigenvalue calculations that are exact in the chain length for oper- 
ators with more than two i?-charge impurities. Furthermore, since the goal is to eventually 
compare anomalous dimensions with 1/J energy corrections to corresponding string states 
near the pp-wave limit of AdSc, x S 5 , and because the string angular momentum J is pro- 
portional to the lattice length L, this virial expansion is precisely what is needed to devise a 
practical method for testing the AdS / CFT correspondence at any order in the gauge theory 
loop expansion for an arbitrary number of i?-charge (or worldsheet) impurities. 

The purpose of this paper is to provide a detailed description of these methods, and 
to compare previously derived numerical predictions ^U] with corresponding predictions 
from Bethe ansatz technology near the BMN limit. We will focus on three particular closed 
sectors of the theory, each labeled by the subalgebra of the full superconformal algebra which 
characterizes the spin variables of the equivalent spin chain system. Specifically, there are 
two sectors spanned by bosonic operators and labeled by su(2) and sl(2) subalgebras plus 
an su(2|3) sector which includes fermionic operators. Section 2 is dedicated to an analysis of 
the bosonic su(2) closed sector to three-loop order in A. In section 3 we analyze an 5u(l|l) 
subsector of the closed su(2|3) system to three- loop order. The spin-chain Hamiltonian in 
the bosonic st(2) sector has previously been determined to one loop, and we analyze this 
system in section 4. We conclude in the final section with a discussion of future applications 
and directions of study. 

2 The 5u(2) sector 

Single-trace operators in the closed su(2) sector are constructed from two complex scalar 
fields of N = 4 SYM, typically denoted by Z and 0. Under the SO (6) ~ U(l) R x SO (A) 
decomposition of the full S77(4) i?-symmetry group, the Z fields are charged under the 
scalar U{1)r component and is a particular scalar field carrying zero i?-charge. The basis 
of length- L operators in the planar limit is constructed from single-trace monomials with / 
impurities and total .R-charge equal to L — I: 

Tr^Z 1 - 1 ) , Tr(0 J - 1 Z0Z z '- / - 1 ) , Tr^^Z^Z 1 - 1 ' 1 ) , ... (2.1) 

The statement that this sector of operators is "closed" means simply that the anomalous 
dimension operator can be diagonalized on this basis, at least to leading order in large N c 

The heart of the spin-chain approach is the proposition that there exists a one- dimensional 
spin system whose Hamiltonian can be identified with the large- N c limit of the anomalous 
dimension operator acting on this closed subspace of operators ^2]- Since the anomalous 
dimensions are perturbative in the 't Hooft coupling A, it is natural to expand the su(2) spin 
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chain Hamiltonian in powers of A as well: 

n ^ ' 

Comparison with the gauge theory has shown that successive terms in the expansion of the 

(2) 

Hamiltonian have a remarkably simple structure: the one- loop-order Hamiltonian H s i 2 s is 
built out of permutations of pairs of nearest-neighbor fields and, at n-th order, the Hamil- 
tonian permutes among themselves fields which are at most n lattice sites apart. This is a 
universal structure which leads to remarkable simplifications in the various closed sectors of 
the theory [T3J. 

Beisert, Kristjansen and Staudacher [3] have introduced the following useful notation for 
products of permutations acting on operators separated by an arbitrary number of lattice 
sites: 

L 

{ni,n 2) ...} = ^Pk+ 

ni,fc+ni+l-^fc+ri2, £.'+71.2+1 " > 

(2.3) 

fc=i 

where Pjj simply exchanges fields on the i th and j th lattice sites on the chain. The upshot of 
the gauge theory analysis is that the equivalent spin-chain Hamiltonian for the su(2) sector 
can be written in a rather compact form in terms of this notation. The result, correct to 
three-loop order, is (see jS] for details) 

H iu(2) = 2({}-{0» (2.4) 
H %) = 2(-4{} + 6{0}-({0,l} + {l,0})) (2.5) 
H i%) = 4 ( 15 " 26 {°} + 6 ({°» !} + {!> 0» + {°> 2 } " ({°- 1, 2} + {2, 1, 0})) . (2.6) 

The form of the three-loop term H^l^ was first conjectured in [5. based on integrability 
restrictions and BMN scaling; this conjecture was later corroborated by direct field-theoretic 
methods in |7j (see also for relevant discussion on this point). Our goal is to develop 
practical methods for finding the eigenvalue spectrum of the spin-chain Hamiltonian for 
various interesting cases. 

2.1 One- loop order 

We start at one-loop order with H^ 2 s in eqn. (j2.4j) . which provides a natural 'position- 
space' prescription for constructing matrix elements in an J-impurity basis of operators. As 
an explicit example, we consider first the basis of two-impurity operators of length L = 8: 

Tr(0 2 Z 6 ) Ti((f)Z(f)Z 5 ) Tr((f)Z 2 (j)Z 4 ) Tr(0Z 3 0Z 3 ) . (2.7) 
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It is easy to see that the one-loop Hamiltonian mixes the four elements of this basis according 
to the matrix 



U su(2) 
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2 -2 








\ 




-2 4 


-2 









-2 


4 






\ 





-2^2 
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/ 



(2- 



This tri-diagonal matrix generalizes to arbitrary L and it is simple to show that the two- 

r(2) 
( su(2) 



(2) 

impurity one-loop eigenvalues of H L, are given by the formula [To] 



^(2) _o•2/ 7^r^ ^ „_ n „ _/(^- 2 )/ 2 > ^ even 



^. = 8 -%r3iJ -»•••••«— iU-sK i odd- w> 

Although we defer our discussion of the Bethe ansatz approach until later in this section, 
we note that the two-impurity su(2) Bethe equations for this spin chain can be solved 
exactly and their solution agrees with eqn. ()2.9|) . Note that if the denominator L — 1 were 
replaced by L, the above expression would agree with the usual lattice Laplacian energy for 
a lattice of length L. The difference amounts to corrections to the free Laplacian of higher 
order in 1/L and we will seek to understand the physical origin of such corrections in what 
follows. 

To compare gauge theory predictions with 1/ J corrections to the three- impurity spectrum 
of the string theory on AdS^ x S 5 , we need to determine the large- L behavior of the three- 
impurity spin chain spectrum. We are primarily interested in systems with few impurities 
compared to the length of the spin chain and we expect that impurity interaction terms in 
the Hamiltonian will be suppressed by powers of the impurity density (i.e. inverse powers 
of the lattice length). This suggests that we develop a virial expansion of the spin chain 
Hamiltonian in which the leading-order term in 1/L gives the energy of free pseudoparticle 
states on the lattice (labeled by lattice momentum mode numbers as in the two-impurity 
spectrum eqn. (|2.9)0 and higher 1/L corrections come from TV-body interactions described 
by vertices Vn- A reasonable guess about how the iV-body interactions should scale with 
1/L suggests that we can write the one-loop-order energy for / impurities in the form 

A 1 n 7T 21 A 
E({ni}) =/+ — ^sin 2 ^ + ^ z ^- T ^ body (n 1 ,...,n / ) + --- , (2.10) 

i=l N=2 

where the leading-order contribution / measures the naive dimension minus -R-charge, the 
next term is the lattice Laplacian energy of I non-interacting pseudoparticles and the 1/L 
corrections account for interactions between pseudoparticles (which may depend on the lat- 
tice momenta mode numbers n$). In the many-body approach, one would try to derive such 
energy expressions by rewriting the Hamiltonian in terms of creation/annihilation operators 
b nv £>t . for the pseudoparticles (commuting or anticommuting as appropriate). The iV-body 
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interaction vertex would generically be written in terms of the b, W as 

N N 

v N = J2 5 

niH \-n N ,mi-\ \-rn N 

ni,m,i i=l i=l 

where /W({^i}, {^«}) is some function of the lattice momenta and the Kronecker delta 
enforces lattice momentum conservation. One has to determine the functions /jv by matching 
the many-body form of the Hamiltonian to exact spin chain expressions such as eqn. ()2.4|) . 
We will see that, once the Hamiltonian is in many-body form, it is straightforward to obtain 
a density expansion of the higher-impurity energy eigenvalues. 

The discussion so far has been in the context of one-loop gauge theory physics, but the 
logic of the virial expansion should be applicable to the general case. To include higher-loop 
order physics we must do two things: a) generalize the functions /jv({^i}, { m i}) defining the 
multi-particle interaction vertices to power series in A and b) allow the free pseudoparticle 
kinetic energies themselves to become power series in A. We will be able to carry out the 
detailed construction of the higher-loop virial Hamiltonian in a few well-chosen cases. To 
match this expansion at n-loop order in A to the corresponding loop order (in the modified 
't Hooft coupling A' = g\ M N c / J 2 ) in the string theory, we need to determine the Hamiltonian 
to 0(L~( 2n+1 )) in this virial expansion. (The first curvature correction to the pp-wave string 
theory at one loop, for example, appears at 0(X'/ J) or, in terms of gauge theory parameters, 
at 0(X/L 3 ).) Auspiciously, it will turn out that this virial expansion in the su(2) sector is 
truncated at small orders in 1/L, allowing for simple eigenvalue calculations that are exact 
in L (although perturbative in A). 

The first step toward obtaining the desired virial expansion is to recast the spin chain 
Hamiltonian H su ^), which is initially expressed in terms of permutation operators, in terms 
of a creation and annihilation operator algebra. We begin by introducing the spin operators 

S ± = l -{a x ±ia y ) S z = ^a z , (2.12) 

where a are the Pauli matrices and S^, Sj act on a two-dimensional spinor space at the j th 
lattice site in the chain. In this setting the Z and fields are understood to be modeled by 
up and down spins on the lattice. The nearest-neighbor permutation operator P^i+i can be 
written in terms of spin operators as 

Pi,i+i — S^S i+1 + S i + 2S*S? +1 + - , (2.13) 
and the one-loop Hamiltonian in eqn. ()2.4j) can be written as 

L L 

H<£ m = ~ £ (S;SJ +1 + Sr S ^ - 2 S-SUi + ~ 2 ■ (2-14) 

j=l 3=1 

A Jordan- Wigner transformation can now be used to express the spin generators in terms 
of anti-commuting creation and annihilation operators (anti-commuting because each site 
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can be either unoccupied (Z) or occupied once (0)). A pedagogical introduction to this 
technique can be found in [TH] - The explicit transformation is 

S+ = b]K(j) = K(J)b) 
SJ = K(j)bj b,Kij) 

S? = b]b s - 1/2, (2.15) 



where the Klein factors 



3-1 

K{j)=exp[iTr^2blb k ] (2.16) 



k=l 



serve to ensure that spin operators on different sites commute, despite the anticommuting 
nature of the bj. The functions K(j) are real, Abelian and, for j < k, 

[K(j),S k ]=0. (2.17) 

The operators bj and bj can therefore be written as 

b) = S+K(j) bj = SjK{j) , (2.18) 

and we easily verify that they satisfy the standard anticommutation relations 

{bj, b\} = 8 jk {b],bl} = {b 3 ,b k } = 0. (2.19) 

Cyclicity on the lattice requires that S^+i = Si, a condition which can be enforced by the 
following boundary condition on the creation and annihilation operators: 

L 

b L+1 = (-l) m 6i I = J2b f jbj, (2.20) 

3=1 

where the integer / counts the number of spin chain impurities. Since we are primarily 
interested in the three-impurity problem (for comparison with the corresponding string re- 
sults reported in JU]), we will henceforth impose the boundary conditions in eqn. (|2.20j) for 
odd impurity number only. We can use all of this to re-express eqn. (J2.14JI in creation and 
annihilation operator language, with the result 

L 

H %) = E + h U h ** - b U b i - b h+i + 2 fy] + iWs+i) ■ ( 2 - 21 ) 

3=1 

Converting to momentum space via the usual Fourier transform 

6 3 - = ^=ge^6 p (2.22) 



sfL „ 

v p=0 
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yields 



L-l 



H 



(2) 
flu(2) 



4 E 



sm 



p=0 



L-l 



+ 7. £ 



27ri(g- 



(2.23) 



p,g,r,s=0 



This is a rather standard many-body Hamiltonian: it acts on a Fock space of momentum 
eigenstate pseudoparticles, contains a one-body pseudoparticle kinetic energy term and a 
two-body pseudoparticle interaction (the latter having the critical property that it conserves 
the number of pseudoparticles). Note that the Hamiltonian terminates at two-body interac- 
tions, a fact which will simplify the virial expansion of the energy spectrum. This termination 
is a consequence of the fact that the one-loop Hamiltonian contains only nearest-neighbor 
interactions and that lattice sites can only be once occupied. 

Because the pseudoparticle (or impurity) number is conserved by the interaction, three- 
impurity eigenstates of the Hamiltonian must lie in the space spanned by 



ki + ^2 + ks = mod L 



(2.24) 



where the ground state \L) is identified with the zero-impurity operator Tr(Z L ) and the 
condition of vanishing net lattice momentum arises from translation invariance on the spin 
chain (which in turn arises from the cyclicity of the single-trace operators in the operator 
basis). As a concrete example, the basis of three- impurity states of the L = 6 su(2) spin 
chain is 

mil \l) bim \l) mi\ \d vm i^> , (2.25) 

and the matrix elements of the Hamiltonian (J2.23)) in this basis are easily computed: 

\ 



H 



(2) 
SU(2) 



V 



1 



1 


1 


3 


3 


-1 


-1 


19 


1 


3 


3 


1 


19 


3 


3 



(2.26) 



The first-order perturbation theory corrections to the three-impurity operator anomalous 
dimensions are the eigenvalues of this matrix. 

The construction and diagonalization of the Hamiltonian matrix on the degenerate basis 
of three-impurity operators can easily be carried out for larger L. The results of doing this 1 
for lattice sizes out to L = 40 are displayed in figure [T] According to eqn. 1)2.10)) . we expect 

(2) 

the eigenvalues of H^L, to scale for large L according to 



^E^\{k}) + ±E^({h}) + 0(\L^) 



(2.27) 



1 Using the position- or momentum-space formalism is purely a matter of convenience. In practice we 
have found that for all sectors the momentum-space treatment is computationally much more efficient. The 
large-L extrapolations of both methods can be checked against each other, and we of course find that they 
are in agreement. 
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10 20 30 40 

L 

Figure 1: One-loop su(2) spin chain spectrum vs. lattice length L (6 < L < 40) 

The scaling coefficients E^?X an< l ^W2) can eas ily be extracted from the data displayed in 
figure H by fitting the spectral curves to large-order polynomials in 1/L (a similar treatment 
was used in [H]). The results of this procedure are recorded for several low-lying levels in 
the spectrum (excluding zero eigenvalues) in table [T] As originally reported in [TU] , string 
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£(1,3) /£(1,2) 


Lattice Momenta (k\, ^2, ^3) 


1 + 2.6 x 10" 


-9 


2 


-4.9 x 10- 
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2 - 5.0 x lO" 7 


(1,0,-1) 


3 + 4.6 x 10- 


-9 
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-8.8 x 10- 
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7/3-3.0 x 10~ 7 


(1,1,-2) 


3 + 4.6 x 10- 


-9 


7 


-8.8 x 10- 


7 


7/3-3.0 x 10~ 7 


(-1,-1,2) 


4 + 6.0 x 10- 


-9 


8- 


- 1.1 x 10- 


6 


2 - 2.9 x lO" 7 


(2,0,-2) 


7 + 3.2 x 10- 


-8 


14 


- 7.1 x 10 


-6 


2 - 1.0 x 10~ 6 


(1,2,-3) 


7 + 3.2 x 10- 


-8 


14 


- 7.1 x 10 


-6 


2 - 1.0 x 10~ 6 


(-1,-2,3) 


9 + 2.2 x 10- 


-7 


18 


- 5.1 x 10 


-5 


2 - 5.7 x 10" 6 


(3,0,-3) 


12 + 5.7 x 10 


-5 


28 


+ 3.8 x 10 


-3 


7/3 - 1.4 x 10~ 3 


(2,2,-4) 


12 + 5.7 x 10 


-5 


28 


+ 3.8 x 10 


-3 


7/3 - 1.4 x 10~ 3 


(-2,-2,4) 


13-5.6 x 10 


-5 


26 


-3.8 x 10 


-3 


2 + 1.3 x lO- 3 


(1,3,-4) 


13-5.6 x 10 


-5 


26 


-3.8 x 10 


-3 


2 + 1.3 x 10" 3 


(-1,-3,4) 



Table 1: Scaling limit of three-impurity su(2) numerical spectrum at one loop in A 
theory makes the following simple predictions for the large-L su(2) expansion coefficients 




E 



.(1.2) 
su(2) 



{k\ + k\ + kl)/2 

„(1,3) , „(1,2) _ 2 
- C/ su(2)/- C/ su(2) — Z 

p (l,3) /p (l,2) _ 7 
^(2)/ ^(2) — 3 



fci + fc 2 + h = 



(*i ^ h ^ k 3 ) 



(fa = k 2 , h 



2fci) • 



(2.28) 



Note the slight annoyance that we must distinguish the case where all mode indices are 
unequal from the case where two indices are equal and different from the third. The last 
column of tabled displays the choice of indices {ki} that best fit each spectral series and the 
other columns display the deviation of the extrapolation coefficients from the string theory 
predictions of eqn. IJ2.28]) . As the lattice momenta increase, higher-order 1/L corrections to 
the spectrum become stronger and more data will be required to maintain a given level of 
precision of the polynomial fit. This effect can be seen directly in the extrapolated eigenvalues 
in table ^ Nonetheless, it is clear from the table that the gauge theory match to the string 
theory prediction is extremely good. 

We also note that the spectrum in table [T] exhibits a degeneracy of eigenstates whose 
momentum labels are related by an overall sign flip (a symmetry that is implemented on 
the operator basis by a parity operator P which reverses the ordering of all fields within the 
trace) . This degeneracy among "parity pairs" of gauge theory operators was observed in [5| , 
where it was shown that it arises as a consequence of integrability (which can, in turn, be 
used to constrain the form of the Hamiltonian at higher loop order |Hj ) . See (2H| for further 
discussion on the implications of this degeneracy. 

To corroborate these results on the gauge theory side we turn to the one-loop Bethe 
ansatz for the Heisenberg spin chain. The Bethe ansatz for chains of spins in arbitrary rep- 
resentations of arbitrary simple Lie groups was developed some time ago [21] (see also [22] 
for an extension to supersymmetric spin chains) and applied only recently to the specific case 
of the dilatation operator of M = 4 SYM [T2*l IT5] . In the notation of ^3], the Bethe equa- 
tions are expressed in terms of the so-called Bethe roots (or rapidities) U{ associated with 
the various impurity insertions in the single-trace ground state TrZ L . In a one-dimensional 
dynamical interpretation, the impurities are pseudoparticle excitations and the roots param- 
eterize in some fashion the lattice momenta of the pseudoparticles. The index i in the Bethe 
root Ui runs over the total number / of impurities. A second index qi = 1, . . . , 7 is used to 
associate each of the I Bethe roots with a particular simple root of the sl(4|4) symmetry 
algebra associated with M = 4 SYM. The Bethe ansatz then takes the form (see JT3j and 
references therein for further details) 



where V qi denotes the q* Dynkin coefficient of the spin representation and M is the Cartan 
matrix of the algebra. To be slightly more specific, if a qi are the root vectors associated with 




L 



i 



(2.29) 
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the nodes of the Dynkin diagram and /x is the highest weight of the spin representation, then 
the Dynkin coefficient (for a bosonic algebra) is V q . = 2a' ft ' • \i / '(a^) 2 and the elements of 
the Cartan matrix are M quq . = 2 ■ / (a^) 2 (note that diagonal elements M quqi = 2). 
(For superalgebras see, eg., [2HJI2I]) Furthermore, since the spin chain systems of interest 
to us are cyclic and carry no net momentum (analogous to the level-matching condition in 
the string theory), the Bethe roots Ui are subject to the additional constraint 



n 



(2.30) 



Finally, having found a set of Bethe roots Ui that solve the above equations, the corresponding 
energy eigenvalue (up to an overall additive constant; see, eg., [13J) is given by 



E = T, -T^TTT • (2.31) 




In the current application all impurities are of the same type (i.e. carry the same Dynkin 
label), so the index qi can be ignored. It is worth noting, however, that the Dynkin coefficient 
V qi can vanish, in which case the associated Bethe roots do not contribute directly to the 
energy. 

The Bethe equations are typically exactly soluble for the case of two identical impurities 
(i.e. two Bethe roots Ui, ui associated with the same simple root of the algebra). The 
two-impurity su(2) Bethe equations, for example, yield solutions that reproduce the familiar 
two-impurity anomalous dimension formula noted above in eqn. ()2.9j) (see ^21 IS] for further 
examples). For three and higher impurities, however, exact solutions are not known. Since 
we are ultimately interested in comparing with string theory predictions at large values 
of the S 5 angular momentum J, an alternate approach is to solve the Bethe equations 
perturbatively in small 1/L. Experience shows that, in the limit where we can neglect 
interactions between excitations (or impurities), the Bethe roots are simply the inverse of 
the conserved momentum carried by the impurities. With a little work, one can show that 
the Bethe ansatz conditions, eqns. (J2.29I2.30|) . can be solved order-by-order in a large- L 
expansion: 

u t = (L + A i VZ + B i + ---) , (2.32) 

where < ki < L is the usual integer lattice momentum. The half-integer powers of L may or 
may not be present in eqn. ()2.32|) : they are needed to deal with special kinematic situations 
(such as when a pair of impurities have the same lattice momentum) where the integral 
power expansion would be singular. The eigenvalues of the spin chain (or the anomalous 
dimensions of the corresponding gauge theory operator) are then obtained as a power series 
in 1/L by substituting the expansion of the Bethe roots into eqn. (j2.31|) . This is the approach 
introduced by Minahan and Zarembo for the so (6) spin chain in [T2*| . Since we wish to carry 
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out similar calculations at higher orders in A, we will review this methodology at one-loop 
order for the specific case of three identical impurities in the su(2) spin chain. (Since the 
su(2) chain is a subsector of the so (6) system studied in ^2], the three-impurity Bethe ansatz 
predictions derived here are of course implied by the all- impurity so (6) anomalous dimension 
formula derived in jT2| at one loop.) 

We now apply this to the closed su(2) sector where the Dynkin diagram has a single node, 
the Cartan matrix is M su ( 2 ) = 2 and the Dynkin coefficient of the fundamental representation 
is Vsu(2) = 1- Consequently, the Bethe equations ()2.29I2.30|) reduce to 

u i + i/2\ L 'fu 1 -u 1±1 \ 

Ai. \ ?/• — it ■ — i I 



Ui — i/2j *-?-\Ui — Uj 



n 



Ui + i/2 
Ui — i/2 



(2.34) 



With three or more pseudoparticle excitations, bound-state solutions can arise which satisfy 
the second equation (|2.34j) . These solutions are characterized as having pseudoparticle states 
sharing the same lattice momenta (eg. fcj = kj for the i th and j th roots). The generic solutions 
to the Bethe equations can therefore be loosely divided into those which do or do not contain 
bound states. For three impurities with no bound states present (k% ^ k 2 ^ k 3 ), eqn. ()2.34|) 
states that k% = —k\ — k 2 . The strategy of [T^j can then be used to obtain a systematic 
expansion of su(2) Bethe roots in powers of L , with the result 

L-A 3A4 ntr- 1 ) 

Ul ~ 2nki + ix(ki - k 2 ){2k 1 + k 2 ) + } 

{L-A)k\ + {L-A)k l k 2 -2{L-l)kl 
U2 2nk 2 {k 2 1 + k l k 2 -2k 2 2 ) + 1 ' 

(L-l)kl-(8-5L)k 1 k 2 + 2(L-l)kl . 

2n(k 1 + k 2 )(2k 1 + k 2 )(k 1 + 2k 2 ) +U{L >■ [Z - 6b) 

Substituting these roots into the energy formula eqn. (|2.3U|) gives the following expression 
for the anomalous dimension of the su(2) three-impurity operator at one-loop: 



2 

E %)( k ^ k 2) = ^{k\ + hk 2 + kl) {L + 2) + 0(L" 4 ) (h + k 2 ^ k 3 ) . (2.36) 



This is in perfect agreement with the string theory results of eqn. f!2.28|) and, of course, the 
numerical gauge theory results in table ^ When a single bound state is present the Bethe 
roots must be altered. Taking, for example, k\ = k 2 , the cyclic constraint in eqn. ()2.34|) sets 
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&3 = —2ki, and the Bethe roots are 

= -7 + 3n/L + 3L 

u 2 = _ 7 + 3^-3L + ; 
D7r/ei 

In this case the anomalous dimension is predicted to be 

E %)( k i) = JT k K 3L + 7 ) + °(^ 4 ) (*i = fc 3 = -2^) , (2.38) 

which is again in agreement with the results of eqn. ()2.28|) and tabled (note that the fractional 
powers of L~ l have obligingly canceled out of the final expression for the energy). 



2.2 Two and three-loop order 

A similar analysis can be performed on the two-loop su(2) spin-chain Hamiltonian. As before, 
we use the Jordan- Wigner transformation restricted to an odd-impurity basis of operators 
to rewrite the two-loop Hamiltonian (J2.5j) in terms of position-space fermionic oscillators, 
obtaining a result similar to eqn. (J2.21j) : 

l , , 

-3 6}6 i -46}6} +1 6 i 6 i+ i 

j=l y 

+b) +l b) +2 b 3 b 3+1 + b]b] +1 b j+1 b j+2 + &J&t +! M- +2 J . (2.39) 
Passing to momentum space, we obtain the two- loop analogue of eqn. (|2.23j) : 



H 



(4) 
SU(2) 



J j+2 



b, + b% 



+2 



bj + b]b 3 . 



h$» = -sE- 4 (t)^ 

i L ~ l 

1 V— \ / 2m(q + r) -2-Ki(p+s) 4iTi(q-s) 27ri(g-s) \ ~ 

+Y \ e L +e L + e L -4e £ J b\b\b r b s 8 p+q>r+s . 

p,q,r,s=0 

(2.40) 



Although the two-loop Hamiltonian includes "long-range" interactions among non-neighboring 
lattice sites, the momentum-space Hamiltonian ()2.40|) conveniently terminates at two-body 
interaction terms. An equally important point is that, for fixed momenta p,q, . . ., the one- 
body (two-body) operators scale as L~ A {L~ 5 ) for large L (the corresponding scalings for the 
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one-loop Hamiltonian were L~ 2 (L~ 3 )). This special relation between density scaling and 
power of coupling constant is critical for matching to string theory. 

We deal with the problem of finding the eigenvalues of the combined one- and two- 
loop Hamiltonian via Rayleigh-Schrodinger perturbation theory: at each value of the lattice 
length L we treat the one-loop operator H^L^ as a zeroth-order Hamiltonian and regard 

fl^ 2 j as a first-order perturbation. The 0(X 2 ) corrections to the spectrum of fl2 2 j are then 

found by taking expectation values of the perturbation H^ 2 ^ in the (numerically-determined) 

(2) 

eigenvectors of H^Ly This is the recipe for non-degenerate first-order perturbation theory 
and we might worry that the previously-noted parity-pair degeneracy of the eigenvalues of 
^iu(2) wou ld force us to use the rules of degenerate perturbation theory. As discussed in 
nHll2n]> however, parity degeneracy can be traced to the existence of a higher Abelian charge 
which is conserved to at least three-loop order. This charge can be used to show that the 
formulas of non-degenerate perturbation theory can be used without modification. The basic 
observation is that conservation of the Abelian charge guarantees that the matrix element of 
^su(2) between two degenerate eigenstates of H^L* with different eigenvalues of the higher 
Abelian charge vanishes: this eliminates the vanishing energy-denominator singularities that 
would otherwise invalidate the non-degenerate first-order perturbation theory formulas (and 
similar arguments apply to the higher-order cases). 

Using this method, we have evaluated the 0(A 2 ) corrections to the spectrum of anomalous 
dimensions for lattice sizes from L = 6 to L = 40. As before, we fit the spectral data to a 
power series in 1/L to read off the leading scaling coefficients of the low- lying eigenvalues. 
As mentioned in the discussion of the two-loop Hamiltonian ()2.40j) . we expect the two- loop 
eigenvalues to have the following scaling behavior in 1/L: 

E?{{h}) = ^E^({h}) + ^E^({h}) + 0(A 2 L- 6 ) . (2.41) 

The numerical data confirm that the eigenvalues scale at least as fast as L -4 . The resulting 
numerical values for the leading scaling coefficients of low- lying eigenvalues, E^fl and E^S., 
are presented in table 121 As originally reported in JU] > string theory makes the following 
simple predictions for the two-loop large- L expansion coefficients: 

E f u §) = ~(k 2 + k\ + klf/lQ h + k 2 + k 3 = 



E^/E^ = y (h = k 2 , h = ~2h) . (2.42) 

The low-lying levels in the table match the string theory predictions quite accurately and 
the decline in precision as one goes to higher energies is expected. As a consistency check we 
note that this time we have no freedom to choose the momenta (ki, k%, ks) associated with 
each state: they have been fixed in the one-loop matching exercise. 

The three-loop su(2) Hamiltonian ()2.6|) can be dealt with in a similar fashion. The posi- 
tion space operator version of this Hamiltonian is too long to record here, but its momentum 
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£(2,4) 




£(2,5) 
su(2) 




7+2,5) / 7+2,4) 
^(2)/ rj su(2) 


(k 1 ,k 2 ,k 3 ) 


-0.25-4.6 x 10~ 9 


-2 + 8.0 x 10- 


7 


8 - 3.4 x 10~ 6 


(1,0, 


-1) 


-2.25 - 1.4 x 1(T 6 


-19 


+ 2.6 x 10 


-4 


76/9 + 1.2 x 10~ 4 


(1,1, 


-2) 


-2.25 - 1.4 x 1(T 6 


-19 


+ 2.6 x 10 


-4 


76/9 + 1.2 x 10" 4 


(-1,- 


-1,2) 


-4 + 8.3 x 1(T 7 


-32 


- 1.1 x 10 


-4 


8 + 3.0 x 10- 5 


(2,0, 


-2) 


-12.25 - 9.9 x 1(T 6 


-98 


+ 2.3 x 10 


-3 


8 - 2.0 x 10" 4 


(1,2, 


-3) 


-12.25 - 9.9 x 1(T 6 


-98 


+ 2.3 x 10 


-3 


8 - 2.0 x 10" 4 


(-1,- 


-2,3) 


-20.25 + 3.2 x 10~ 3 




-161.4 




7.97 


(3,0, 


-3) 


-36 - 2.8 x 10~ 3 




-304.6 




8.46 


(2,2, 


-4) 


-36 - 2.8 x 10~ 3 




-304.6 




8.46 


(-2,- 


-2,4) 


-42.25 + 4.9 x 10~ 3 




-337.0 




7.97 


(1,3, 


-4) 


-42.25 + 4.9 x 10~ 3 




-337.0 




7.97 


(-1,- 


-3,4) 



Table 2: Scaling limit of three-impurity su(2) numerical spectrum at two loops in A 
space version is fairly compact: 



H 



(6) 
su(2) 



L-l 1 L-l , 

Ef. /VH\ ~l~ 1 ii 27vi(q+r) 27Ti(2q+r) 27T»( 9 +2r) 27vi(q-3s) 

sin 6 lt-Ulb p + — 2^ < -10 + + + 



p=0 



p,q,r,s=0 



2iri(q — r — 2s) 

10 e E 



2iri(2q — 2r-3s) 27ri(3q-2r — 3s) 2_7ri(£— r— 3s) 2Tri(2q-r — 3s) 2TTi(q-2s) 

+e l + e l + e l + e l _ e z 

2ni(2q — r — 2s) 2ivi(3q — r — 2s) 27ri(g+r — 2s) 2ivi(q — s) 47vi(q — s) 67vi(q — s) 

—e l _ e i _ e e + 29 e 1 — 10 e 1 + e E 



27Ti(2g-s) 27Tt(3q-3) 27ri(q+i — s) 27ri(2q+r-s) 27rt(q + 2r-s) ~, ~,~ ~ 

-e L +e L -e ^ +e £ +e ^ }blblb r b x 5 n 



J p u q u r u s Up+q jr + s 



+ 



1 L_1 
f2 Yl 



p,q,r,s,t,u=0 



2iri(q+3r-2t-3u) 2ifi(q + 2r- s-2t-3u) 

e l + e £ 



27Ti(2 9 +3r-t-3n) 2^(q+2r + s-tt) , ~ ~ 

fe z + e z ^J&J&JfcAM 



(2.43) 



It contains at most three-body operators and a careful examination of terms shows that, for 
fixed momenta, the one-body operators scale as i+ 6 , the two-body operators as L~ 7 and 
so on. We therefore expect the leading scaling coefficients in the 0(A 3 ) eigenvalues to be 



(3,6) 
«u(2) 



and E^flj to use a by-now-familiar notation. To find the eigenvalues to this order, we 



E 

continue with the Rayleigh-Schrodinger perturbation theory strategy: the 0(A 3 ) correction 
to any eigenvalue is the sum of the matrix element of H^L, in the appropriate eigenvector 
of h[ 2 }^ plus the second-order sum-over-states contribution of 



su(2) j jau -° ua - l ^ o^^f iiu-ui uci ouni-u v ci-oiai^o uuiiiiiuuiiun ui ±J ~ s\x(2}' These two pieces 
can easily be computed numerically from the explicit Hamiltonian operators at a fixed L. 
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Parity degeneracy and conservation of the higher Abelian charge mentioned above continue 
to hold, and we can again use non-degenerate perturbation theory formulas to compute the 
eigenvalue corrections. We have generated numerical eigenvalue data for lattices from L = 6 
to L = 40 and the large- L scaling coefficients of the low-lying states extracted from those 
data are given in table El As is by now well-known, the detailed match to string theory 





^u(2) 


p (3,7) /p (3,6) 
• C/ «j(2)/' C/ ai(2) 


1 7„ Ir. 

[ki, fc s 


, fa) 


0.1250 


2.0003 


16.003 


(1,0, 


"I) 


4.125 


58.03 


14.07 


(1,1, 


-2) 


4.125 


58.03 


14.07 


(-1,- 


•1,2) 


7.999 


128.2 


16.03 


(2,0, 


-2) 


49.62 


713.3 


14.37 


(1,2, 


-3) 


49.62 


713.3 


14.37 


(-1,- 


-2,3) 


91.15 


1,454 


15.96 


(3,0, 


-3) 


263.8 


3,739 


14.17 


(2,2, 


-4) 


263.8 


3,739 


14.17 


(-2,- 


•2,4) 



Table 3: Scaling limit of three-impurity su(2) numerical spectrum at three loops in A 

breaks down at three-loop order, so there is no point in trying to match these results to 
string predictions. We were initially motivated to pursue a virial treatment of this problem 
because the more classic Bethe ansatz methods were not yet able to deal with non-local spin 
chains. Important progress has recently been made on the Bethe ansatz side and it may be 
useful to compare what can be done by each method (results should of course agree). 

A modified Bethe ansatz for the su(2) sector of the gauge theory, possibly incorporating 
all orders of higher-loop physics, has recently been proposed jSJ E] 2 It is an instructive 
exercise, and a useful consistency check on this bold proposal, to verify that it reproduces 
the higher-loop scaling coefficients for three impurity anomalous dimensions that we have 
just computed by virial methods (and displayed in tables [TJ |2]and|IH). F° r completeness, we 
briefly summarize the new ansatz, referring the reader to jHj for a detailed account. In the 
new ansatz, the momenta pi of the excitations (closely related to the Bethe roots) become 
functions of A (as well as L and mode numbers) and are determined by a modified version 
of eons. (12.3312.341) : 

= T\ ^\-^\ +l V> = 0. (2.44) 

Dependence on A enters through the phase function ip(pi) which is defined in terms of the 

2 The long-range ansatz based on the Inozemtsev spin chain in 8 suffers from improper BMN scaling at 
four- loop order, a problem that is surmounted in [5]. 
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excitation momenta pi as follows: 



1 , /. A 



cpfa) = ^ cot (pi/2) y 1 + ^ sin 2 (ft/2) . (2.45) 

The energy eigenvalue corresponding to a particular root of these equations is given in terms 
of the excitation momenta pi by the formula 




£ s „, 2 1 = ^ — ! \ / 1 + ^ sin 2 ( ft /2) - 1 ) ( 2 . 16) 



Finding exact solutions of these equations is even more difficult than before, but we can follow 
the previous strategy of developing an expansion in powers of 1/L about non-interacting 
impurities on an infinite lattice. This is achieved by expanding the excitation momenta pi 
according to 

L n=l L 2 

where the integers ki (subject to the cyclicity constraint Yli h = 0) characterize the non- 
interacting state about which the expansion is developed. The appearance of half-integer 
powers of L _1 in this expansion is needed to accommodate bound-state solutions to the 
Bethe equations which arise when some of the momenta ki are equal. Solutions to the Bethe 
equation 1)2.44)1 will determine the expansion coefficients p\ in terms of the mode numbers 
ki and ultimately lead to expansions of the energies as power series in L -1 , with coefficients 
that are functions of X/L 2 . 

Explicit results for the L _1 expansion of gauge theory operators of arbitrary impurity 
number, derived by the above method, were presented in |T3]. 3 As usual, expressions are 
different depending on whether all momenta are unequal or some subset of them are equal. 
For all mode numbers ki unequal the J-impurity energy formula in [T5\ is 



E su (2) 




(2.48) 



where we have used A' = A/ J 2 = A/(L — I) 2 for convenience (J = L — I is the total R- 
charge). To compare with our virial results, we must further expand in A; expanding to first 

3 It is important to note that the focus of this paper is a different Bethe ansatz, designed to match the 
spectrum of the string theory: the gauge theory Bethe ansatz results are derived for comparison purposes. 
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and second order yields the following scaling coefficients (valid for all ki unequal): 

E M2) = k\ + hk 2 + k\ = 2(k\ + hk 2 + kl) 

= +QT + r 2 ) 2 E™ = -2(q 2 + qr + r 2 ) 2 . (2.49) 

These one- and two-loop coefficients match the numerical results presented in tables Q and 
Eland the analytic string formulas in eqns. (J2.28p2.42j) . It is harder to write down a general 
formula for the many cases in which subsets of momenta are equal but the solution for the 
particular case of three impurities with a two-excitation bound state (ki = hi = n, k% = —2n) 
was also presented in [To] : 

E su{2 ) = L - 3 + 2Vl + X'n 2 + Vl + A' An 2 

\'n 2 ( 1 6 12 8 \ , 

+ ^^ + ^T-T7—, -^ ^ , o -( 2 - 5 °) 



L-3V.1 + A'n 2 Vl + X'n 2 Vl + A'4n 2 ^1 + X'n 2 s/l + X'An 2 , 
To compare with the virial results, one must again expand the energy in powers of A. Doing 
so yields the following one- and two-loop bound-state scaling coefficients: 

= -j" 4 *S?, = -is* 4 - (2.5i) 

We easily verify that this agrees with numerical virial results to two-loop order. 

The three-loop coefficients obtained by expanding the energy formulas in eqns. ()2.48I2.50|) 
are given by 

,(3,6) I 



E fn{2) = ( 2 fc i 6 + 6 k i 5 k 2 + 15 h 4 k 2 2 + 20 k x 3 k 2 3 + 15 h 2 k 2 A + 6k x k 2 5 + 2 k 2 6 ) 

(2.52) 



£§2) = \ ( 8 k i + 24 k i k 2 + 51 h A k 2 2 + 62 h 3 k 2 3 + 51 h 2 k 2 A + 24 k x k 2 5 + 8 k 2 6 ) 



for (ki 7^ k 2 7^ fc 3 ), and 



= f« 6 ^«?,=58« 6 , (2.53) 



for the bound-state solution with (ki = k 2 = n, fc 3 = —2n). The numerical values of these 
0(A 3 ) coefficients are tabulated for several low-lying states in the spectrum in tabled The 
correspondence with table El which displays the three-loop expansion coefficients extracted 
from numerical diagonalization of the three-loop Hamiltonian, is good. At this order in the 
loop expansion higher-order l/L corrections to the spectrum are more important (compared 
to the one- and two- loop cases), and the numerical extrapolation is less reliable (especially 
as the lattice momenta increase). The precision can always be improved by including data 
from larger lattices in the extrapolation. We emphasize that this discussion concerns the 
different methods of calculation of operator dimensions in the su(2) sector only. It seems to 
us to give useful further evidence that the long-range Bethe ansatz for the su(2) sector of 
the gauge theory jHj is exact. 
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/o r-\ 

^u(2) 


/O IT S 

E (3,7) 


p (3,7) /p (3,6) 
£/ «u(2)/- C/ «l(2) 


(/ci,/c 2 ,/c 3 ) 


0.125 


2 


16 


(1,0,-1) 


4.125 


58 


14.06 


(1,1,-2) 


4.125 


58 


14.06 


(-1,-1,2) 


8 


128 


16 


(2,0,-2) 


49.625 


713 


14.37 


(1,2,-3) 


49.625 


713 


14.37 


(-1,-2,3) 


91.125 


1,458 


16 


(3,0,-3) 


264 


3,712 


14.06 


(2,2,-4) 


264 


3,712 


14.06 


(-2,-2,4) 



Table 4: Three-impurity su(2) spectrum from the long-range Bethe ansatz at three loops 

3 A closed su(l|l) subsector of su(2|3) 

The three-impurity string theory analysis of ^U] identified a fermionic sector of the theory 
which is diagonalized by string states composed of fermionic excitations projected onto par- 
ticular four-dimensional subspaces (which transform in an SU(2) 2 x SU(2) 2 notation as a 
(2, 1; 2, 1) or (1, 2; 1, 2) of SO(A) x SO (A)) and symmetrized in their SO(A) x SO (4) indices. 
It was also shown that this three-impurity subsector of the theory decouples at all orders in 
A. 

On the gauge theory side this subsector corresponds to an su(l|l) subgroup of the closed 
su(2|3) sector studied by Beisert in pHE]. (Supersymmetric integrable su(n\m) spin chains 
have previously been studied in certain condensed-matter applications; see, eg., pH].) In 
the present setting the fields of su(2|3) consist of three complex scalars (p a and two complex 
fermions ip a . In the closed su(l|l) subspace we restrict to a single scalar denoted by Z and a 
single fermion labeled by if). Just as in the su(2) sector, we use the fermionic position-space 
oscillators bp bj to create or annihilate fermionic if> insertions in a ground state composed 
of L scalars: 

\L) = Tr(Z L ) b) \L) = Tt{Z 1 ■ ■ ■ Z^Z j+x ■ ■ ■ Z L ) . (3.1) 

In j7j, Beisert gave the action of the Hamiltonian on the su(2|3) spin chain to three- 
loop order. 4 In the notation of jjj, the action of the Hamiltonian on basis states can be 
represented in terms of special permutation operators denoted by 

(A,...A N ) 

\B 1 ...B N j ' 

4 Beisert's three-loop Hamiltonian was restricted in [J] to the bosonic sector, but the author has since 
provided us with the complete version. 
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which replace all occurrences of the upper sequence of fields A\ . . . An in the trace by the 
lower sequence B\. . . B N . Restricting Beisert's su(2|3) Hamiltonian to the su(l|l) subsector 
at one-loop order yields 



In terms of the position-space oscillators of eqn. (j3.1|) . the su(l|l) Hamiltonian can be as- 
sembled by inspection and takes the form 

L 

H Sm = E ( b h + - b U b J ~ • (3-3) 

3=1 

There are no higher-body interaction terms at this order in A. This fact can be checked by 
computing 

(L\b l+ MHS m )b\bl +1 \L) = 2 , (3.4) 

which reproduces the two-body matrix element given by the last term in eqn. (|3.2j) . In 
momentum space we obtain 

<V)= 4 X> a (f) s &- (3 - 5) 

p=0 

The two-loop su(l|l) momentum-space Hamiltonian can be extracted in the same manner 
(the position-space version is too long to print here): 

L-i h-\ f 

2-ni(q-2r) 2ni(2q — r) 2iri(q-r) 



p=0 p,q,r,s=0 ^ 



e l + e l — 4e l 



2iri(q — 2r — s) 2-7rz(q+s) 2iri(q — r + s) 2iri{2q — 2r — s) ■ ~ , ~ . — ~ 

-2 e z - 2 e—T— + e z + e z - W&}& r 6 a S p+q , r+s . (3.6) 
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Finally, the complete three-loop Hamiltonian for this subsector is 

L-l . L-l f 

( 6 S ^-^ . R /»7T\ ~ + ~ 1 V"^ 60vri(q-r) 2^(27q-29r) 27vi(2Sq-29r) 

^u(i|i) = 32 E sm (t)^ p "16 S |2e i +2e z 

2Tri(27g-28r) 2jri(29q-28r) 27ri(29g-27r) 567ri(g-r) 587ri(qr-r) 

-4 e z + 37 e e 6 e z + 8 e z — - 72 e z — 

27ri(29<j-29r-2s) 27ri(29q-30r-s) 27ri(29q-29r-s) 27ri(29<3-28i — s) 

-6 e 2 - 40 e z + 37 e e _ 8 e e 

27ri(27g-28r+s) 27ri(28q-28r+s) 27ri(29qr-28r+s) 27ri(27q-27r+s) 

+8 e e + 2 e e - 40 e e - 4 e e 

27ri(29g-27r+s) 27H(27q-27r+2s) 2ni(29q-30r-2s) | ~, ~, ~ ~ 

+8 e e + 2 e e + 8 e e W p b\b r b a 5 p+q , r+s 
1 L - x r 

1 ^ — > 27ri(g + 2i — 3a — 2t) 27ri(g+3i — 3s — 2t) 2iri(q+2r — 3s-t) 

H > <^2e e - e e _4 e e 

16 ^ 1 

p,q,r,s,t,u=0 v 

27rz(2g + 3r-3s — t) 2ivi(q + 2r — 2s — t) 2iri(2q+3r — 2s — t) 2iri(q + 2i — 3s — It — u) 

-e e + 8e e + 2e e - 4 e e 

27ri((j+37 — 3s-2t — u) 2iri(q+2i — 2s+u) 

+2 e e + 2 e e 

27ri(q + 2r— s + n) 27ri(q+2r — 2s - t+u) I ~ ~ 

-4 e e _ 4 e e ^j>T&;& s & t &« 5 p+q+r>s+t+u . (3.7) 

We note that -£^u(i|i)' ^iufili) anc ^ ^iufiii) terminate at one-body, two-body and three-body 
interactions, respectively. This will permit us to obtain the exact L-dependence of successive 
terms in the A expansion of energy eigenvalues. 

As in the su(2) sector, we can use non-degenerate perturbation theory to extract the L^ 1 
scaling coefficients of the su(l|l) eigenvalue spectrum up to three loops in A. The scaling 
coefficients extrapolated from numerical diagonalization of lattices up to L = 40 are recorded 
for one-loop, two-loop and three-loop orders in tables HI and [7J respectively. The same 
increase in leading power of L' 1 with corresponding order in A that was noted in the su(2) 
sector is found here as well (we use the same notation for the scaling coefficients as before 
in order to keep track of these powers). It should also be noted that, because the impurities 
in this sector are fermions symmetrized on all group indices, the lattice momenta of all 
pseudoparticles must be different. The string theory results of [TU] amount to the following 
predictions for the one-loop and two-loop scaling coefficients: 

^fi]i) = (*? + + k\) E^ m = 

< 4 i 3 - \tft + hh + klf = -{k{ + hk 2 + klf . (3.8) 

The agreement of these predictions with the data in tables El and El is excellent (with the 
usual caveat that data on larger and larger lattices is required to maintain a fixed precision 
as one goes to higher and higher energy levels). 

The scaling limit of the three- loop ratio ^wi|i)/-^Wi|i) is recorded for the first few low- 
lying states in the spectrum in table [71 These values are in disagreement with the corre- 
sponding three-loop predictions from the string theory as can be seen by comparing with the 
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^u(l|l) 




^udli) 




p (l,3) /p (l,2) 


{k 1 ,k 2 ,k 3 ) 


1 - 


- 1.3 x 10- 


10 


-1.9 x 10 


-8 


-1.9 x lO" 8 


(1,0, 


-1) 


4- 


- 1.0 x 10- 


-7 


1.8 x 10- 


5 


4.6 x 10~ 6 


(2,0, 


-2) 


7- 


- 2.5 x 10" 


-7 


4.4 x 10- 


5 


6.3 x 10~ 6 


(1,2, 


-3) 


7- 


- 2.5 x 10~ 


-7 


4.4 x 10- 


5 


6.3 x 10~ 6 


(-1,- 


-2,3) 


9- 


-3.9 x 10" 


-7 


7.9 x 10- 


5 


8.7 x 10- 6 


(3,0, 


-3) 


13 


-4.0 x 10 


-6 


8.2 x 10- 


4 


6.3 x 10~ 5 


(1,3, 


-4) 


13 


-4.0 x 10 


-6 


8.2 x 10- 


4 


6.3 x 10~ 5 


(-1,- 


-3,4) 


16 


- 2.0 x 10 


-5 


4.1 x 10" 


3 


2.6 x lO" 4 


(4,0, 


-4) 


19 


-3.5 x 10 


-5 


7.3 x 10- 


3 


3.8 x 10" 4 


(2,3, 


-5) 


19 


-3.5 x 10 


-5 


7.3 x 10- 


3 


3.8 x 10" 4 


(-2,- 


-3,5) 



Table 5: Scaling limit of one- loop numerical spectrum of three- impurity su(l|l) subsector 



£(2,4) 
s«(l|l) 


P (2,5) 


p (2,5) /p (2,4) 
^SU(1|1)/ ^.611(1 11) 


(fci, k 2 , h) 


-0.25 


-0.99999 


3.99995 


(1,0,-1) 


-4.00006 


-15.990 


3.998 


(2,0,-2) 


-12.251 


-48.899 


3.992 


(1,2,-3) 


-12.251 


-48.899 


3.992 


(-1,-2,3) 


-20.25 


-80.89 


3.995 


(3,0,-3) 


-42.25 


-168.2 


3.98 


(1,3,-4) 


-42.25 


-168.2 


3.98 


(-1,-3,4) 


-64.00 


-254.6 


3.98 


(4,0,-4) 


-90.26 


-359.3 


3.98 


(2,3,-5) 


-90.26 


-359.8 


3.99 


(-2,-3,5) 



Table 6: Scaling limit of two-loop numerical spectrum of three-impurity su(l|l) subsector 

results of ^U]. Given the well-established three-loop disagreement between the string and 
gauge theory in the su(2) sector, however, this disagreement in the su(l|l) subsector is not 
unexpected. 

The extrapolated gauge theory results in eqn. (|3.8|) for the one-loop coefficients Kl, 

(1 2) 

and -Egu^jx) should be checked against the predictions of the general one-loop Bethe ansatz 
|12[ ITB*] applied to the su(l|l) sector (as far as we know, no higher-loop Bethe ansatz is 
available here). To apply the general Bethe ansatz equation of eqn. (|2.29|) . we note that the 
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p(3,7) / p (3,6) 
-^11(1 1 1) / ^011(1 11) 


(fci, k 2 , h) 


-86.41 


(1,0,-1) 


-85.71 


(2,0,-2) 


-83.74 


(1,2,-3) 


-83.74 


(-1,-2,3) 


-101.9 


(3,0,-3) 


-96.01 


(1,3,-4) 


-96.01 


(-1,-3,4) 


-158.1 


(4,0,-4) 



Table 7: Scaling limit of three-loop numerical spectrum of three-impurity su(l|l) fermionic 
subsector 

su(l|l) Dynkin diagram is just a single fermionic node: the Cartan matrix is empty and the 
single Dynkin label is Vwiii) = 1 |2*3*1 124j . We therefore obtain the simple one-loop Bethe 
equation 




Rather remarkably, eqn. ()3.9|) can be solved exactly for arbitrary impurity number! The 
general su(l|l) Bethe roots are 

1 [ kiii\ 

Ui = 2 cot \17) ( ^ 

and the energy eigenvalues computed from eqn. (|2.31|) are 

^u(i|i)=4^sin 2 (^) , (3.11) 

with the usual condition ki — mod L from eqn. ()2.30|) . This is just the sum of free lattice 
Laplacian energies and clearly matches the energies one would obtain from the one-loop 
su(l|l) Hamiltonian of eqn. (J3.5|) (since the latter has no interaction terms). No expansion 
in 1/L was necessary in this argument, but it is straightforward to expand the energies in 
1/L and verify the numerical results obtained in table El and eqn. ()3.8|) . 



4 The sl{2) sector 

As noted in JU], integrable s((2) spin chains have previously been the subject of several 
studies involving, among other interesting problems, high-energy scattering amplitudes in 
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non- super symmetric QCD (see, eg., and references therein). The sl(2) closed sector of 
M = 4 SYM was studied in ^7j, and the spin chain Hamiltonian in this sector is presently 
known to one loop in A. 

The constituent fields in this sector are SO (6) bosons Z carrying a single unit of i?-charge 
(Z = 05 + and each lattice site on the sl(2) spin chain is occupied by a single Z field 
acted on by any number of the spacetime covariant derivatives T> = T>\ + iT> 2 . The total 
i?-charge of a particular operator is therefore equal to the lattice length L, and an /-impurity 
operator basis is spanned by single-trace operators carrying all possible distributions of I 
derivatives among the L lattice sites: 

Tr {V : Z Z L - 1 ) , Tr (V^Z VZ Z L ~ 2 ) , Tr (l? 1 ^ ZVZ Z L ~ Z ) , ... (4.1) 

The integer / counts the total number of derivatives in the operator and, since any number 
of impurities can occupy the same lattice site, one can think of n derivative insertions at the 
i th lattice site as n bosonic oscillator excitations at the i th lattice position: 

(a]) n \L) ~ Tr (Z i ~ 1 'D n ZZ L - i ) , ... (4.2) 

The ground state \L) is represented by a length L chain with no derivative insertions: \L) = 
Tr(Z L ). 

The one-loop sl(2) spin chain Hamiltonian (corresponding to the dilatation operator in 
this sector) was constructed in [T7j and was defined by its action on basis states rather than 
directly expressed as an operator: 

L 

rr(2) _ ST^ ttB\{2) 
Sl(2) - 2-^/ 3,3+1 ' 
J = l 



H^\a\y{a^\L) = Y J 

j'=0 



5 j=j/ (h(j) + h{n - j)) 



\j-r\ 



{a\Y{a\)^'\L) 

(4.3) 



(where h(n) = 1 + . . . + 1/n are the harmonic numbers). In other words, H S J 2 ) ls a sum 
over the position-space Hamiltonian //fj+i which acts on the j th and (j + l) th (neighboring) 
lattice sites; the action of Hfj+i can be summarized by the explicit form given for 
above. Since it is only defined by its action on the state {cl\) 3 \a\) n ~ 3 \L), it is difficult to 

(2) 

immediately translate H^J^ to momentum space. However, it is possible to expand it in 
powers of fields and use eqn. ()4.3j) to iteratively determine the expansion coefficients. The 
virial argument furthermore tells us that higher powers in the fields will determine higher 
powers of L _1 in the expansion of the energy. For our current purposes, it suffices to know 
the Hamiltonian expanded out to terms of fourth order in the fields and this truncation of 
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the Hamiltonian can easily be constructed by inspection: 

L 



H 



(2) 

m 



3=1 



2a) + a)_, 



1 t 2 



,t 2 



o t 2 , t 2 

2a) +a)_ l 



+ 



(4.4) 



Transformation to momentum space gives 

L-l 



p=0 



L-l 



4 £ * 



p+q,r+s 



sin 



2^ 



. 2 

sm — 

LJ 



2 (p + g)^ , t _ _ 

• I n\n\n n 

L 



+ sm ) apa' q a r a s + 



(4.5) 



p,q,r,s=0 

This Hamiltonian acts on an J-impurity Fock space spanned by the generic states 



(4.6) 



with lattice momenta labeled by ki = 0, . . . , L — 1, and subject to the constraint ki = 
mod L. Numerically diagonalizing this Hamiltonian on a range of lattice sizes, we obtain 

(1 2) 

data from which we extract the numerical predictions for the one-loop coefficients E^JJ 

E s xJ presented in table |H1 String theory makes the following predictions [Ej for the see 
coefficients 



and 



E 



E 



(1,2) 
0[(2) 

(1,2) 
S [(2) " 



(kl + hk 2 + k\ 



E 



(1,3) / P (l,2) 



3n 2 



p (l,3) /P (1,2) 
js[(2) / ^(2) 



•1(2) 

-7/3 



-2 h^k 2 ^ h 

k 2 = n,k 3 = -2n , 



(4.7) 



and we can easily verify that the agreement with table [H] is excellent. 

The extrapolated predictions can again be checked against those of the corresponding 
one-loop Bethe ansatz equations. In the sl(2) sector the highest weight is —1/2: the Dynkin 
diagram therefore has coefficient V s ^2) = — 1 and the Cartan matrix is M S U2) = 2. The Bethe 
equations ()2.29J2.3()|) thus reduce to 



U; 



i/2 



Ui + i/2 
1 = 



lit 



n 



n 

3& 

Ui — i/2 



Uj + i 



U ; 



(4.8) 
(4.9) 



Ui +i/2 t 

Apart from a crucial minus sign, this is identical to the su(2) Bethe equation ()2.34|) . In the 
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£(1,2) 
^[(2) 


£(1,3) 
^[(2) 




p (l,3) /F (1,2) 
^51(2) / -^[(2) 




1 + 1.2 x 1(T 9 


-2 - 3.1 x 10- 


7 


-2 - 3.1 x 10- 7 


(1,0, 


-1) 


3 - 7.6 x 10- 9 


-7+1.9 x 10- 


6 


-7/3 + 6.3 x 10- 7 


(1,1, 


-2) 


3-7.6 x 1(T 9 


-7 + 1.9 x 10- 


6 


-7/3 + 6.3 x 10- 7 


(-1,- 


-1,2) 


4-2.8 x 10~ 7 


-8 + 6.9 x 10- 


6 


-2 + 1.7 x 10- 6 


(2,0, 


-2) 


7- 2.9 x 10- 7 


-14 + 7.1 x 10- 


-5 


-2 + 1.0 x 10- 5 


(1,2, 


-3) 


7- 2.9 x 10~ 7 


-14 + 7.1 x 10- 


-5 


-2 + 1.0 x 10- 5 


(-1,- 


-2,3) 


9-4.1 x 10" 7 


-18 + 1.0 x 10- 


-4 


-2 + 1.0 x 10- 5 


(3,0, 


-3) 


12 + 8.4 x 1(T 7 


-28 - 1.5 x 10 


-4 


-7/3 - 1.2 x lO" 5 


(2,2, 


-4) 


12 + 8.4 x 10~ 7 


-28 - 1.5 x 10 


-4 


-7/3 - 1.2 x 10~ 5 


(-2,- 


-2,4) 


13-7.0 x 10~ 6 


-26 + 1.7 x 10- 


-3 


-2 + 1.3 x 10- 4 


(1,3, 


-4) 


13-7.0 x 10~ 6 


-26 + 1.7 x 10- 


-3 


-2 + 1.3 x 10- 4 


(-1,- 


-3,4) 


16 - 1.4 x 10~ 6 


-32 + 3.9 x 10- 


-4 


-2 + 2.4 x 10~ 5 


(4,0, 


-4) 


19 - 7.5 x lO" 6 


-38 + 2.2 x 10- 


-3 


-2 + 1.1 x 10" 4 


(2,3, 


-5) 


19 - 7.5 x 10~ 6 


-38 + 2.2 x 10- 


-3 


-2 + 1.1 x 10~ 4 


(-2,- 


-3,5) 


21-3.4 x 10" 6 


-42 + 8.8 x 10- 


-4 


-2 + 4.2 x lO" 5 


(1,4, 


-5) 


21-3.4 x 10" 6 


-42 + 8.8 x 10- 


-4 


-2 + 4.2 x 10- 5 


(-1,- 


-4,5) 



Table 8: Scaling limit of numerical spectrum of three-impurity st(2) sector at one-loop in A 
absence of bound states, eqn. ()4.8|) is satisfied by the following Bethe roots: 

2(1 + L)k\ - (4 + L)k x k 2 - (4 + L)k\ 



2^k x {k\ + k x k 2 - 2k 2 ^ 



0{L~ l ) 



2(1 + L)kj- (4 + L)hk 2 - (4 + L)kj 
2 2Ttk 2 (k\ + k x k 2 -2kl) + 1 ' 

2{l + L)kj + {% + hL)k l k 2 + 2{l + L)kl 
U3 2n(k 1 + k 2 )(2k 1 + k 2 )(k 1 + 2k 2 ) +U{L } - l4 - iUJ 

Using eqn. ()2.3H) . we obtain 

E %)(ki, k 2 ) = A [kl + kl k 2 + kl) (L - 2) + 0(L- 4 ) {h + k 2 ^ k 3 ) . (4.11) 
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For the bound state characterized by k\ = k 2 = n and k 3 = 



— 2n, the Bethe roots are 



7-3^ + 3^ 



+ 0(L 




Ml = 



67m 



7 + 3^ + 31, 



+ 0(L 




u 2 = 



6nn 



u 3 = 




(4.12) 



with spin chain energy 




(3L - 7) + 0(L- 4 ) {h = k 2 = n, k 3 



2n) . 



(4.13) 



These results again agree with the numerical results in table |S] (and match the corresponding 
three-impurity string theory predictions in [TU]). 

5 Conclusions 

We have demonstrated that the virial expansion of the Af = 4 SYM spin chain Hamilto- 
nian for small impurity number provides a simple and reliable method for computing exact 
anomalous dimensions of mult i- impurity operators at small scalar i?-charge (chain length) 
and estimating with great precision the near-BMN scaling behavior of these dimensions as 
the .R-charge becomes large. The latter application, which is suited to direct comparison 
of gauge theory predictions with corresponding results on the string side of the AdS/CFT 
correspondence, works well for three-impurity operators to three-loop order in A in the 
su(2) sector (the order to which the su(2) Hamiltonian is known definitively). Specifically, 
the numerical predictions from the virial approach for the near-BMN scaling coefficients 



( E m$V E M2)> E iSv £ S?)' E SS) and E m&) match corresponding results from the su(2) 



long-range Bethe ansatz to three-loop order, and agree with near-plane-wave string theory 
predictions to two loops (the disagreement with string theory at three loops is by now an 
expected outcome in these studies). We also find convincing agreement near the BMN limit 
between the virial approach and the Bethe ansatz results at one-loop order in the closed 
sl(2) and su(l|l) subsectors. As a side result we have found in the su(l|l) sector an exact (in 
chain length) agreement between the Bethe ansatz and the virial expansion for one-loop op- 
erator dimensions with arbitrary impurity number (this was only possible because the Bethe 
equations can be solved exactly in this subsector for any number of impurities). There are 
currently no higher-loop Bethe ansatze for the sl(2) and su(l|l) systems, however, so in this 
sense our numerical predictions go beyond the current state of Bethe ansatz technology (see 
[2*7] for further developments of higher- loop gauge theory physics in non-su(2) sectors). It 
would be very interesting to find a general long-range Bethe equation appropriate for J\f = 4 
SYM at higher loop-order in A, both for comparison with string predictions and with the 
virial approach studied here. 
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